THEORETICAL AND EMPIRICAL ANALYSES 



OF 

THICK SKINNED SWEPT WING PANELS 



by 



Paul Whitney Parcells 
Lieutenant, United States Navy 
B.S., United States Naval Academy, 1960 



Submitted in partial fulfillment of the 
requirements for the degree of 

MASTER OF SCIENCE IN AERONAUTICAL ENGINEERING 

from the 

NAVAL POSTGRADUATE SCHOOL 
December 1967 



(L- - v C \ 



ABSTRACT 

The validity of conventional skin and stringer type 
analyses is verified experimentally for a thick milled skin 
aircraft structure. The accuracy of solution is evaluated 
by comparing the strain energies in a mathematical model to 
that obtained experimentally from a wing panel of a high 
performance aircraft. 

Results indicate that th| conventional methods of a 
skin-stringer/matrix force approach are valid for a thick 
milled skin aircraft structure of this type provided the 
panel loading is neither severe nor concentrated. Addi- 
tionally the grid elements chosen must be small enough to 
preclude any gross averaging errors. 



/ 



I 



DUDLEY KNOX LIBRARY 
NAVAL POSTGRADUATE SCHOOL 



TABLE OF CONTENTS 



CHAPTER PAGE 

I. INTRODUCTION 13 

II. THEORETICAL ANALYSIS 16 

III. DISCUSSION OF THEORETICAL ANALYSIS 28 

IV. EXPERIMENTAL ANALYSIS 30 

V. DISCUSSION OF EXPERIMENTAL ANALYSIS 46 

VI. COMPARISON OF RESULTS 48 

VII. CONCLUSIONS AND RECOMMENDATIONS 52 



REFERENCES 
APPENDIX A 
APPENDIX B 
APPENDIX C 
APPENDIX D 
APPENDIX E 



54 

THEORETICAL DETERMINATION OF SHEAR FLOWS 56 



EQUIPMENT CALIBRATION CURVES 61 
STRAIN GAGE LOCATIONS 62 
STRAIN GAGE INSTRUMENTATION 64 
PRINCIPAL STRESSES AND ORIENTATIONS 66 



3 



LIST OF TABLES 



TABLE 

I. 

II. 

III. 



Theoretical Panel Parameters 
Theoretical Strain Energy Distribution 
Average Sub-Panel Experimental Stresses 



PAGE 

22 

27 

42 



5 



LIST OF FIGURES 



FIGURE 






PAGE 


1 . 


Sketch of F8U-3 Wing Center Section 




17 


2. 


Center Section of F8U-3 Port Wing 




18 


3. 


Idealized Panel on F8U-3 Port Wing 




20 


4. 


Panel Showing Dimensions for Edge Loading Ratios 


23 


5. 


Exploded View of Idealized Wing 




25 


6. 


Sketch of Strain Gage Locations (Outer 


Skin) 


34 


7. 


Sketch of Strain Gage Locations (Inner 


Skin) 


35 


8. 


Sketch of Planform of Partitioned Panel 


41 


9. 


Skin Thickness Variation with Span 




43 


10. 


Experimental Principal Stresses 




45 



7 



LIST OF PLATES 



PLATE 

1. F8U-3 Wing Test Installation 

2. Wing Mounting for Load Application 

3. Interior Strain Gage Installation 

4. Equipment Installation 



TABLE OF SYMBOLS 



A 

a n 



a . 

i 



a c 

AFT-IB 

b 

CIB 

E 

FB 

FIB 

FWD-IB 

G 

1 b 

1 d 

q 

RB 

RIB 

T 

tj. 



U 

V 

V„ 



Inclosed wing cell area, square inches 
Panel width (tip) , inches 

Panel width (root) , inches 
Panel end loading (tip) , pounds 
Panel end loading (root) , pounds 
Aft intermediate beam 
Panel length, inches 
Center intermediate beam 

Modulus of Elasticity, pounds per square inch 
Front beam 

Front intermediate beam 

Forward intermediate beam 

Shear modulus , pounds per square inch 

Panel span (forward) , inches 

Panel span (aft) , inches 

Shear flow, pounds per inch 
Rear beam 

Rear intermediate beam 

Torque, inch pounds 

Panel thickness (tip) , inches 

Panel thickness (root) , inches 
Internal strain energy, inch pounds 
Constant voltage source, volts 
Wheatstone Bridge output voltage, volts 



11 



w i 

w c 

Xw 

^w 

p 

8 

/*- 

8 

cr 

8 

© 

0p 



Panel edge length, inches 
Panel edge length, inches 
x-axis co-ordinate of wing, inches 
y-axis co-ordinate of wing, inches 

Panel angle (AFT-IB to pivot rib) , degrees 

Panel angle (RIB to intermediate rib) , degrees 

Poisson's ratio 

Strain, micro-inches per inch 

Normal stresses, pounds per square inch 

Shear stresses, pounds per square inch 

Angle of twist per unit length, radians per inch 

Angle of principal axis measured from the first 
gage of the rosette, degrees 



12 



CHAPTER I 



INTRODUCTION 

Aircraft structures in general, and wings in parti- 
cular, have been constructed using the conventional stringer 
and sheet methods. This form of construction, by nature not 
a simple structure, introduced the problem of how to accu- 
rately predict the load distribution and especially the 
stress concentrations that could be developed. The oldest 
representation scheme for the discrete element analysis of 
the problem evolved as the axial force member-shear panel 
idealization that is described in Reference 1. This method 
is particularly acceptable to the aircraft industry, as it 
readily lends itself to various matrix techniques that are 
easily adaptable to a digital computer for solution. Gen- 
erally speaking, this method of solution makes no provision 
for accounting for Poisson's ratio and sweep effects. 

With the advent of the very heavy milled skins of 
modern high performance aircraft, one immediately questions 
the validity of the conventional form of solution. In order 
to properly account for bending-torsion interaction, chord- 
wise curvatures, sweep coupling, shear deflections, and the 
increasing use of large cutouts, several alternative methods 
of analysis have been developed. 

The most widely used solution technique is outlined by 
J. H. Argyris and S. Kelsey in Reference 2. The aircraft 
structure is molded into a large assemblage of elastic ele- 
ments, which permits matrix formulation of the solution in 
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terms of the different energy theorems. The fundamental 
principles originated by Maxwell and Mohr of virtual forces 
and virtual displacements are used in the energy method of 
solution. A practical application of this technique is 
included in References 3 and 4. 

The objective of this study is to determine the valid- 
ity of this lumped parameter form of a theoretical analysis 
for thick skinned structures. Because of the general energy 
approach to the solution, the method chosen to determine the 
accuracy of the results was to compare the total internal 
strain energy of an idealized panel with that of an actual 
thick milled skin aircraft structure of the same dimensions. 
The structure would have to be thick enough to be able to 
support any bending-torsion interaction that may exist in 
addition to the inplane loadings. 

The objective of this report was to be achieved by 
comparing the theoretical strain energy of a structure with 
that determined empirically from strain gage readings. The 
structural specimen to be analyzed was a single wing panel 
taken from the port wing upper skin of a Mach 2.5 all- 
weather fighter project which was canceled by the U. S. 

Navy before the structural test program was completed. 

The necessary strain gage rosettes and accompanying 
instrumentation were installed, and a comparison of the 
theoretical and empirical analyses was made. 

A literature search did not reveal any reference that 
made a comparison of the energy level in an idealized 
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structure to that of an actual wing panel element. The 
importance of a comparison of this type is apparent when one 
realizes that the integrated behavior of a structure is 
dependent upon the energy accounting in the system. 

This study was completed during the 1967-1968 academic 
year at the Naval Postgraduate 8chool, Monterey, California. 
Acknowledgement is gratefully made to Professor C. H. Kahr 
of the Aeronautics Department for his guidance and consulta- 
tion as thesis advisor; R. A. Besel and T. B. Dunton, 
laboratory technicians; and CAPT D. B. Messerschmidt , USMC 
who, while working on a similar thesis project, aided in 
this study. 
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CHAPTER II 



THEORETICAL ANALYSIS 

The structure chosen for study was the wing center 
section of an F8U-3 type aircraft. The particular panel 
analyzed was the port section bounded longitudinally by the 
pivot and intermediate ribs and laterally by the aft- 
intermediate and rear-intermediate beams. The wing center 
section and specific panel are graphically depicted in 
Figure 1. 

The typical industrial approach to analyzing a struc- 
ture of this nature would be to determine what valid sirm-- 
plifying assumptions could be made in order to create a 
mathematically simple idealized structure. The resulting 
model would consist of axial load carrying bara and plane 
shear carrying constant thickness panels. The idealized 
model was developed and is thoroughly discussed in 
References 3 and 4. 

The load application was accomplished by applying a 
torsional moment at the wing fold rib, about the load 
reference axis (CIB) . Considering the port wing inboard of 
the intermediate rib, this couple must be distributed, 
using elementary torque theory, to the equivalent loads on 
a streamwise cross-section (intermediate rib) at station 
y w = 81.98 (Figure 2). The loading was then divided among 
the six cells of the wing. The shear flow for each cell 
was determined in accordance with Reference 1 by equating 
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SKETCH OF F8U-3 WING CENTER SECTION 



FIG. 2 

CENTER SECTION - F8U-3, PORT WING 
(Main structural elements) 
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Wing station, 




the angle of twist of one cell with the remaining five and 
equating torsional moments as follows: 
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T =f 2, i 



( 1 ) 



The results and calculations of the shear flows were 
calculated in Reference 3 and are included in this report 
as Appendix A with inappropriate sections deleted. Using 
the planform dimensions of Figure 2, the shear flow for each 
cell was determined as oriented perpendicular to the CIB. 

Figure 3 shows the idealized panel section of the wing 
with pertinent dimensions obtained from Reference 5. The 
internal strain energy of this skewed shear panel was 
calculated in accordance with S. J. Garvey (Reference 6) as 
follows : 

"Since by similar triangles, the perpendicular distance 
of any point in a trapezoidal panel to a parallel edge is 
proportional to the depth of the panel at that point 
measured parallel with a side, the well-known result follows 
that the "shear flow" in a trapezoidal panel is inversely 
proportional to the square of the depth. 

To find the strain energy, using these geometric sim- 
plifications of a trapezoid, one must solve the following 
equation: 




( 2 ) 
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FIG. 3 

IDEALIZED PANEL ON F8U-3, PORT WING 



where: 




(3) 



Table I lists the appropriate dimensions and constants 
needed for the theoretical analysis. The volume of the 
plate was determined for the idealized panel from: 



This approach of utilizing the average thickness will ob- 
viously limit the accuracy of the resulting theoretical 
strain energy; however, the objective of this part of the 
analysis is to obtain the strain energy by the same conven- 
tional methods currently being used by industry. 

As previously discussed, the plate edge loadings were 
determined from a linear interpolation using the planform 
dimensions. These shear panel edge loadings were distri- 
buted around the perimeter of the panel, in accordance with 
Reference 4, by first assigning a unit shear force to the 
outer edge of each panel. Considering a typical panel as 
shown in Figure 4 and taking moments about appropriate 
points to ensure the maintenance of equilibrium, the follow- 
ing ratios of edge forces were derived: 



V 




(4) 




(5) 
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TABLE I 



PANEL PARAMETERS 



Parameter 


Value 


Units 


Re ference 


/* 


0.32 




7079-T6 Al 


G 


3.9 x 10 6 


psi. 


7079-T6 Al 
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10.3 x 10 6 


psi. 


7079-T6 Al 


a l 


12.132 


in. 


Ref. 5 


a 2 


15.247 


in. 


Ref. 5 
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56.885 


in. 


Ref. 5 


fc r 


0.62 


in. 


Ref. 5 




0.33 


in. 


Ref. 5 
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0.477 


in. 




a.vG 
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61° 13' 


deg. 


Ref. 5 
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58° 52' 


deg. 


Ref. 5 
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TYPICAL PANEL SHOWING DIMENSIONS FOR EDGE LOADING RATIOS. 



( 6 ) 
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( 7 ) 



The Matrix Force Method is described in Reference 11 
and was used to determine the end edge loadings on the 
panels. These loads were obtained from previous calcula- 
tions in Reference 3. The end edge load for the particular 
panel studied in this report (Figure 5) was then scaled to 
the edge utilized in the strain energy calculation of 
Equation (2) . This scaled load was then divided by the 
length of the same edge. The result was the desired edge 
shear flow (q^ = 325.976 pounds per inch). 

Substitution of the appropriate values from Table I and 
the computed shear flow into Equations (2) , (3) , and (4) 

will yield the total shear panel internal strain energy. 

The resulting value was 33.280 inch pounds. 

The idealized model, being a pure shear panel, was in- 
capable of supporting any axial loadings. For this reason 
the effective cross sectional area of the plate was divided 
evenly and distributed to the adjacent flanges (Figures 3 
and 5) . The axial loads were then assumed to be applied 
directly to the flanges. 

In accordance with Reference 9, sets of additional 
axial forces (dashed forces in Figures 3 and 5) were applied 
to the ends of the flanges. The forces are necessary to 
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account for the additional strain energy stored by the axial 
components of shear flows caused by the sweep effect of the 
non-rectangular panels. It is important to realize that 
these sweep coupling loads are internal forces and would not 
be included in the equilibrium equations of the structure. 

From an energy standpoint, these dashed forces account 
for the additional strain energy stored in the panel. The 
energy is accounted for in this manner because the longi- 
tudinal flanges contain the effective area of the panel. 

The contribution of a uniform bar under linearly vary- 
ing axial forces to the total internal strain energy is 
given in Reference 9 as, 

ii= _!_ 

V £AE 

where P^ and Pj are the axial forces applied to the ends 
of the bar. Integration of Equation (8) yields: 

U=-dt(P?+ RB) 

The flange loadings were determined as a percentage of the 
total load based on the ratio of the effective flange area 
of the panel to the total flange area. 

Substitution of appropriate panel parameters and pre- 
viously determined flange loadings into Equation (9) yields 
the additional strain energies required to complete the 
theoretical energy accounting. The distribution of the 
energies is listed in Table II. 
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TABLE II 



THEORETICAL STRAIN ENERGY DISTRIBUTION 
OF THE IDEALIZED MODEL 


Shear Panel 


33.280 


inch pounds 


Forward Flange 


0.313 


inch pounds 


Rear Flange 


1.067 


inch pounds 


Rib Flange 


0.046 


inch pounds 


TOTAL ENERGY 


34.706 


inch pounds 
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CHAPTER III 



DISCUSSION <$F THEORETICAL RESULTS 

The concept of the idealization, as discussed in 
Reference 11, is simply to assign all the direct stress 
carrying capacity of the structural material to the axial 
force members and the shear carrying capacity to the 
shear panel. 

The idealization yields a theoretical solution to the 
problem that utilizes several different averaging processes; 
however, the objective of this report is to analyze the 
accepted conventional methods currently in use in the 
aircraft industry. 

The distribution of the panel area to the adjacent 

flanges is complicated by the apparent double use of the 

area near the tip of the panel. Several wing analyses have 

described alternative methods of including the rib flange in 

12 

the idealized structure. It is of interest to note that 
the theoretical result of the strain energy determined for 
the rib flange is insignificant when compared to the other 
components of the strain energy in the structure. Because 
of the small energy of the rib flange portion of the plate 
contribution, one could conclude that this flange and simi- 
lar pieces do not contribute significantly to the over-all 
wing behavior. Therefore, as discussed in Reference 12, 
analyses of wings with applied normal loadings could be 
simplified considerably by assuming zero flexibility for 
this flange. It should be noted that the application of 
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chordwise loads would require that these elements be in- 
cluded. This approach would appreciably decrease the number 
of redundants in a wing analysis. 

An appropriate study might be the determination of the 
contribution of these rib flanges to the integrated analysis 
of the wing. 
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CHAPTER IV 



EXPERIMENTAL ANALYSIS 

In order to determine the actual strain energy and 
compare it with the value predicted from theory, an experi- 
mental investigation was completed on a panel of an F8U-3 
wing acquired from the Ling-Temco-Vought Corporation. The 
panel section was a typical thick milled skin portion of 
the wing and was tapered along the span as well as in the 
chordwise direction. 

The detailed explanation of the wing installation and 
loading is included in Reference 3; therefore, only the 
highlights of installation and any modifications incorpora- 
ted by the author will be included in this report. 

The wing was mounted inverted on a rigid support jig. 
It was fastened to the jig at the two fuselage pivot lugs 
and the two points on the main box front beam directly 
opposite the bumper points (Plate 1) . Aluminum pads were 
installed under the jig to distribute the loading to the 
laboratory floor. 

The entire wing was subjected to two pure torque loads 
applied at each wing fold rib about the CIB (Pl&te 2) . The 
position of the elastic axis was obtained from manufac- 
turer's data and is shown in Figure 2. The CIB is approxi- 
mately parallel to the elastic axis in the outboard section 
of the wing. It reflects the characteristic breakdown in 
the concept of a load axis toward the root of a swept wing. 
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The predetermined loading (8,000 pounds or 336,000 inch 
pounds of torque) was measured with 10,000 pound dynamo- 
meters connected in series with each hydraulic loading 
cylinder. They were graduated in 100 pound increments but 
could be read accurately to ten pounds. The final load 
magnitudes ware determined with hydraulic manifold pressure 
readings related to pressure-load calibration curves for 
the cylinders found in Appendix B. The load was well 
within the elastic limit of the wing structure but still 
of sufficient magnitude to yield adequate strain levels 
throughout the structure. 

All outer surface strain gage installations on the test 
panel were installed by the author. The internal gages 
were installed by laboratory technicians. All strains were 
measured with SR-4 strain gage rosettes of the type AR-7-2 
manufactured by the Baldwin-Lima-Hamilton Corporation. 
Appendix C is a list of all 102 gages, their gage factors, 
resistances and co-ordinate locations on the panel. The 
rosettes on the outer surface were installed back to back 
with the inner gages in order to give an indication of any 
differential bending effects through the thickness of the 
panel that may have existed. Figures 6 and 7 graphically 
depict the strain gage locations on the panel, and Plate 3 
shows the interior gage installation. A twenty channel 
Budd Company switching and balancing unit was connected to 
an externally powered Wheatstone bridge circuit. The bridge 
and all associated circuitry were an integral part of a 
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Outboard-4 
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FIG. 6 

SKETCH OF STRAIN GAGE LOCATIONS 

(Outer-skin) 



Forward 
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FIG. 7 

SKETCH OF STRAIN GAGE LOCATIONS 
(Inner - skin) 
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Budd/Datran Digital Strain Indicator (Model TC22) that read 
out strain directly in micro-inches per inch. Calibration 
of all equipment preceded each run. The electronic set-up 
is described in Reference 3 and is reproduced in this report 
as Appendix D. 

The procedures used in conducting all test runs were 
as follows: 

1. The digital indicator was calibrated before 
each run ensuring that sufficient time was 
taken between runs to preclude any errors 
being induced from hysteresis effects. 

2. The switching and balancing unit leads were 
plugged into twenty gage terminals at the 
junction panel (Plate 4). 

3. All strain gages and dynamometers were zeroed. 

4. The hydraulic torque loading was applied and 
adjusted to agree with all previous runs. 

5. Equal couple loadings applied to each wing 
tip were ensured by comparing the dynamometer 
readings . 

6. The strain readings were taken and recorded. 

7. The loads were removed. Hysteresis effects 
were removed by allowing sufficient time (8 to 
10 minutes) before commencing the next run. 

Because of ambient laboratory temperature changes, the 
sun striking the wing, and outside electrical interference, 
all test runs were conducted at night. To ensure steady 
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state conditions an adequate (approximately one half hour) 
warm-up time was used for all associated equipment. Repeat- 
ability of strain gage readings was ensured by check runs 
being made periodically and compared with those previously 
obtained. The percentage difference in readings was minute 
and was considered to have an insignificant effect on the 
analysis . 

As stated earlier, the objective of the experimental 
phase of this analysis was to correlate the theoretical 
results that were previously obtained; therefore, once the 
recorded strain data were obtained, the experimental total 
strain energy in the plate was to be determined. Refer *-... -jo 
ence 7 discusses the calculation of strain energy in a 
structure in a general form that is applicable to the panel 
considered in this report. The general form of the total 
strain energy is given by, 



U = fjf [tf + «*)- t ( 5^3 + 



( 10 ) 



+ M + 



0 



Z<=\ 



In the wing panel of this report, where 01 , Tu* , 
and are all nearly zero for plates/panels of this 

type, the stress condition indicated by Equation (10) 
reduces to the following: 

^ < 1J -> 
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This strain energy formula of the stress state in 
Cartesian co-ordinates is the same as that related by 
A. L. M. Grzedzielski in Reference 8. Here the terms have 
the following meanings: Txij term represents the shear 

energy of the panel. The and CT ^ terms reflect 

the inplane and bending/torsion energies. The term contain- 
ing Gy introduces the Poisson's ratio coupling. 

Strain gage data were programmed into an IBM System/360 
Model 67 digital computer to obtain the principal stresses 
and directions at all rosette locations. A FORTRAN IV 
program "RECROSE" was written to convert the strain readings 
to the desired stresses. The discussion of the computer 
program and the results are included in this report as 
Appendix E. 

The test panel was then partitioned into fifteen sub- 
panel elements. The partitioning of the panel (Figure 8) 
would permit a summation of strain energies over the panel 
that would yield the total strain energy. This process is 
a close approximation to the integration required in 
Equation (11) . The strain readings were averaged where more 
than one gage existed. Both the inside and outside rosettes 
on each sub-panel were included. The results of the sub- 
panel averaging are listed in Table III. The sub-panel 
thicknesses used to compute the volumes were the average 
sub-panel thicknesses obtained from Reference 10. A graph 
of the panel thickness versus y w is included in this 
report as Figure 9. 
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FIG. 8 

SKETCH OF PLANFORM OF PARTITIONED PANEL 



TABLE III 



AVERAGE SUB-PANEL STRESSES 


Panel No. 






Tis 


1 


76.94 


-115.93 


-490.25 


2 


119.37 


-173.96 


-338.28 


3 


141.58 


-203.99 


-502.31 


4 


157.43 


-71.62 


-661.05 


5 


237.09 


-104.48 


-669.08 


6 


215.11 


-90.22 


-691. 19 


7 


188.63 


-40.42 


-898.64 


8 


165.23 


-63.82 


-910.20 


9 


244.65 


-80 . 84 


-870.01 


10 


153.65 


-83.44 


-1139.25 


11 


220.78 


-72.57 


-1127.20 


12 


218.42 


-46.80 


-1064.91 


13 


250.56 


- 78.95 


-1225.65 


14 


174.92 


-182.72 


-1291.96 


15 


231. 18 


-106.37 


-1285.93 
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The differences between the inner and outer surface 
stresses (Appendix E) indicate that a bending/torsion effect 
definitely does exist; however, assuming a linear variation 
of the stress through the thickness of the plate produces a 
strain energy that is less than one per cent different from 
a direct stress averaging technique. Therefore, the mid- 
plane stresses obtained by averaging the inner and outer 
surface stresses were used in the calculation of the experi- 
mental strain energies. 

The individual sub-panel strain energies were computed 
from Equation (11) using the required stresses listed in 
Appendix E. The sub-panel energies were then summed to 
obtain the total strain energy of the entire wing panel. 

The result of this summation was 33.714 inch pounds. 

The internal strain energy thus obtained was nearly 
that determined theoretically in Chapter II and was con- 
sidered a reasonably accurate evaluation of the actual 
strain energy that existed in the swept panel element. 

Figure 10 is a plot of the principal stresses obtained 
from "RECROSE" versus location on the panel. The purpose 
of this graph is to give the reader a feeling for the over- 
all panel loading under the applied loading of the whole 
wing. 
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EXPERIMENTAL PRINCIPAL STRESSES 
(Averaged along streamwise chord) 




(rsd) sssjts 
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CHAPTER V 



DISCUSSION OF EXPERIMENTAL RESULTS 

The strain gages on the outer surface were mounted back 
to back with the interior rosettes. The purpose of this 
installation was to detect any bending present that would 
be reflected in a difference in inner and outer strain 
readings. Examination of Figure 10 and Table E4 in 
Appendix E indicates that bending definitely does exist. 
However, when a linear variation of stress through the 
thickness of the plate is assumed, the strain energy does 
not vary appreciably from that computed by direct averaging; 
therefore, even though the wing panel does exhibit the 
effects of bending and twisting, the internal strain energy 
computations do not seem to reflect this effect. The 
experimental analysis averages the inner and outer values 
of stresses before determining the panel strain energy; 
therefore, this method utilizes the midplane stresses that 
are the mean of the two surface stresses. It is significant 
to note that the stresses along the rosette legs must be 
utilized for this computation. The principal stresses 
cannot be used because the orientation of the inner and^ 
outer stresses will not be parallel. 

The experimental result obtained utilizes different 
degrees of an averaging process. Obviously, this will not 
yield precisely the strain energy that actually exists in 
the panel; however, the partitioned panel strain energy 
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that was computed by summing the sub-panel energies must be 
as close as possible to the actual internal strain energy 
that exists. The further partitioning of the panel into 
smaller elements would require the installation of more 
strain gage rosettes. This would possibly further improve 
upon the accuracy of the experimental strain energy results; 
however, this approach was precluded by both time and 
monetary considerations. 

The optimum installation of the strain gage rosettes 
would have been to place them in the center of each of the 
fifteen sub-panel elements; however, the internal gages 
having been previously mounted and the desirability of 
mounting the inner and outer gages back to back precluded 
the optimum installation. Additionally, because of the 
lack of any stress concentrations, the stresses did not 
vary radically over the panel. For these reasons the 
rosette installations were considered adequate for this 
analysis . 
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CHAPTER VI 



COMPARISON OF RESULTS 

The comparison of the theoretical and experimental 
results, obviously, in the final analysis depends on how 
well the energy of the mathematical model approximates the 
true energy of the system. As stated in Reference 8, one 
should not expect a correct result from either a mathema- 
tical model or a lumped averaging process similar to the 
empirical approach taken in this report unless the load 
distribution is uniform over the panel. The strain energy 
error becomes smaller if a more finely partitioned grid is 
used. Thus, when the panel is subdivided into smaller 
elements, the lumped integral form of computing strain 
energy will increase in accuracy. 

The error occurring in the determination of the 
theoretical strain energy is obtained from the following: 

u e " u t 

— = percentage error (12) 

u e 

Utilizing the results obtained in these analyses, the 
percentage error in the theoretical strain energy computa- 
tion was found to be 2.94 per cent. 

It appears that there have been no gross errors made 
in the energy accounting procedures. The comparison made 
indicates a remarkably close agreement between experimental 
and theoretical results. 
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Superficially, one might expect the experimental strain 
energy to be greater than that computed from the idealized 
model. This would seem reasonable when one considers that 
the theoretical solution for a shear panel cannot possibly 
include the strain energies induced from bending, torsion, 
and any cross coupling effects that may exist in the actual 
thick milled skin panel; however, careful scrutiny of the 
idealized model reveals that the Poisson's effect is not 
included; and therefore, there is no mechanism to include 
the chordwise stresses. In fact, they are neglected. On 
the other hand, in the empirical analysis it is shown that 
appreciable chordwise stresses do exist. This is a contra- 
diction to the close agreement (2.94 per cent) in the 
comparison of the strain energies. 

From this discussion one must conclude that the 
stresses determined theoretically in Reference 3 must have 
been in excess of their true values. This excess of span- 
wise and shear stresses would compensate for the lack of 
chordwise stresses. This is the only logical explanation 
for the close comparison of the two analyses. 

In order to determine the validity of this reasoning, 
the principal stresses at station y w = 34.50 inches were 

compared. The stresses in the model were: CT” =714 psi 

m3,x 

and = 596 psi. Similarly, the experimental stresses 

obtained from Figure 10 were: O’"" = 525 psi and 

^ max * 

f max = 517 P si - 
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As was predicted, the idealized normal stresses were 
considerably higher (26.30 per cent) than those obtained 
empirically. Similarly, the theoretical shear stresses were 
higher (13.25 per cent). In general, these overestimates of 
stress values compensated for the lack of chordwise stresses 
when the theoretical strain energy was computed. 

The location of the test panel on the wing did not lend 
itself to the direct application of severe bending and/or 
torsional loads. One can conclude, therefore, from the 
results of this report that when even thick skinned aircraft 
structures are analyzed using the conventional lumped 
stringer and shear panel method, the analysis will be valid 
provided there are no severe bending or torsion loads 
applied directly to the panel. In general, the conventional 
method is valid, provided that the structure is sufficiently 
partitioned into a fine enough grid structure and the load 
distribution is not too irregular. 

In conclusion, one should not expect from an analysis 
based on the lumped stringer and shear panel concept a 
satisfactory agreement with experiment for isolated load- 
ings, particularly when they are represented as point forces 
in a theoretical analysis. In the panel studied and in <■ ' 
other similar cases there will always be ample structural 
material to further partition the model. This partitioning 
of a structure permits the desirable attenuation of the 
effect of load concentrations. These load concentrations 
can create discrepancies between computational and experi- 
mental results. 
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In an associated study, the wing center section is 
being included in the over -all wing analysis. The inclu- 
sion of the energy contribution from the center section 
would not be expected to change the theoretical analysis in 
the vicinity of the particular panel considered in this 
report. This is because of the massive construction of the 
wing in this specific area. Even though the pivot rib is 
not infinitely rigid in the chordwise direction, it is 
nearly so and may be assumed to have zero flexibility. 
Reference 3 verifies this assumption by the close agreement 
between the theoretical and experimental results in the 
vicinity of the particular panel considered. 
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CHAPTER VII 



CONCLUSIONS AND RECOMMENDATIONS 



Conclusions 

It can therefore be concluded that a valid comparison 
was achieved between theoretical and experimental strain 
energies of a thick milled skin aircraft structure. The 
accuracy of the analytical method is subject to the absence 
of concentrated loadings. Additionally, for a theoretical 
structure, a fine element partitioning procedure is 
required. 

The analytical method provided remarkably accurate 
results (2.94 per cent error) and is a valid approach to 
aircraft structures of this specific type. 

Recommendations 

It is recommended that even though the results of this 
study are favorable, another analysis should be made on a 
panel that is subjected to more severe non-planar loadings. 
Possibly a panel outboard of the intermediate rib could be 
used on the same wing. The load would then be applied to a 
thinner section of the wing where there is less torsional 
rigidity, and the point of load application would be closer 
to the panel to be studied. Another possibility would be 
to obtain some means of applying single point loads to the 
same panel at the intersection of the intermediate rib and 
the spanwise beams. Still another method would be to apply 
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the same torque used in this analysis at the intermediate 

rih. 

An additional recommendation, that was previously dis- 
cussed in Chapter III, is that a study be conducted to 
determine the validity of a wing analysis that neglects the 
presence of any rib flanges. 
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APPENDIX A 



THEORETICAL DETERMINATION OF SHEAR FLOWS 



AT SECTIONS PERPENDICULAR TO CIB AT 



y w = 98.7 and y w = 74.3 (Reference 3) 



The theoretical shear flows at sections perpendicular 

to the CIB at y w = 98.70 and y w = 74.30 were calculated in 

order to determine the shear flows at both sections which 

would be used to produce the equivalent applied forces at 

the streamwise intermediate rib, y = 81.98. The cross - 

w 

sections are shown schematically in Figure Al. The thick- 
nesses were averaged mid-panel values, and the web heights 
were taken as the mid-panel height between the upper and 
lower skins. 

A box structure with several cells will have one less 
redundant than the number of cells. In this case, there are 
five redundant webs. It is desired to write six equations 
in six unknown shear flows. This was done by equating the 
angle of twist of one cell with the remaining five, which 
gives five equations, and then writing an equilibrium of 
torsional moments equation. 

The first five equations were obtained by equating the 
angle of twist per unit length, 0 , of one cell with the 

other five using the well-known expression for a box beam, 




(Al) 
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FIG. Al 

SCHEMATIC DIAGRAM OF SECTIONS 
PERPENDICULAR TO THE C.I.B. 
AT y w = 98.7 AND 74.3 





y w =74.3 

Note: dimensions in 
inches 
(ref. 10) 
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The equilibrium of torsional moments may be written as 



T = Z 2 An ^ 



(A2) 



In using Equation A2 the summation is carried out 
around the entire perimeter of each cell. The unit length, 
L, and the constant factor 2G drop out, leaving for cells a 
and b, 




A at t> Abt 



(A3) 



The value of q for any exterior web of cell a is q a ^ 
and for the interior web is (q ~ q^) * Us i n 9 the 
following abbreviations, 



£>kk 

Sak 




AS 



t 




As 

t 



l 




] 



CL— b 



Equation A3 is rewritten, 




(A4) 



The terms &&a. and ^bb represent summations around the 

entire perimeter of their respective cells and & Ja the 
value for the interior web. The enclosed areas were taken 
as the average web height times the distance between webs. 
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Equating the angle of twist per unit length of cell a 



to the remaining cells yields the following equations: 



Making appropriate substitutions, performing the indicated 
algebraic operations and rearranging, one obtains five 
equations in terms of the six shear flows with constant 
coefficients. The required sixth equation was derived as 
previously described. 

The solution to these six simultaneous equations was 
obtained from a digital computer using a FORTRAN program 
which is included in Reference 3. It utilizes Gauss's 
method of elimination with row pivoting and back substitu- 
tion and is designed to yield solutions for one or more 
column vectors forming the right hand side of the set of 
equations. This means that it would produce solutions for 
one or more values of applied torque. The check values of 
q were hand calculated using an iterative procedure 
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suggested by Bruhn in Reference 9 and within slide rule 
accuracy agree well with the computer solutions. 
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Dynomometer reading lb. xIO Error- % of load 



APPENDIX B 

CALIBRATION CURVES FOR DILLON DYNAMOMETERS 
AND HYDRAULIC PRESSURE SYSTEM 
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Hydraulic pressure reading - p.s.i. 




APPENDIX C 



CO-ORDINATE LOCATION OF STRAIN GAGES 
F8U-3 WING PANEL 



Note: Co-ordinate origins are as follows: 

x w = 0 = Intersection of center section droop 
leading edge and the center line of 
the aircraft. 

y TJ = 0 = The aircraft center line, 
w 



UPPER SKIN INSIDE GAGES 



Gage Number 


Gage Type 


X w 


^w 


83-84-85 


AR-7 -2 


137.72 


28.53 


86-87-88 


AR-7 -2 


140.00 


32.25 


89-90-91 


AR-7 -2 


142.55 


30.48 


92-93-94 


AR-7 -2 


145.25 


28.74 


95t96-97 


AR-7 -2 


147.99 


45.40 


98-99-100 


AR-7 -2 


150.30 


43.84 


101-102-103 


AR-7 -2 


152.73 


42.25 


104-105-106 


AR-7 -2 


155.95 


58.68 


107-108-109 


AR-7 -2 


158. 11 


57.21 


110-111-112 


AR-7 -2 


160.24 


55.74 


113-114-115 


AR-7 -2 


163.84 


71.68 


116-117-118 


AR-7 -2 


165.77 


70.51 


119-120-121 


AR-7 -2 


167.72 


69.19 


122-123-124 


AR-7 -2 


169.23 


78.83 


125-126-127 


AR-7 -2 


170.60 


78.83 


128-129-130 


AR-7 -2 


173.00 


78.83 


40-41-42 


AR-7 -2 


168.56 


75. 10 



NOTES 

1. All strain gage rosettes were; mounted at the Naval 
Postgraduate School and oriented with one leg perpendicular 
to the CIB, the 45 degree gage pointing outboard and aft, 
and the third leg parallel to the CIB. 

2. All gages were manufactured by the Baldwin-Lima- 
Hamilton Corporation, taken from lot B-31 and rated at 
120.5 t, 0.5 ohms with a gage factor of 1.97 ± 2 per cent. 

3. The low numbered rosette leg identifies the gage 
and is the leg parallel to the Load Reference Axis. 
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4. Gage No. 98 yields questionable results. 



UPPER SKIN OUTSIDE GAGES 



Gage Number 


Gage Type 


*w 




580-581-582 


AR-7 -2 


137.72 


28.53 


583-584-585 


AR-7-2 


140.00 


32.25 


475-476-477 


AR-7 -2 


142.55 


30.48 


586-587-588 


AR-7 -2 


145.25 


28.74 


589-590-591 


AR-7 -2 


147 . 99 


45.40 


592-593-594 


AR-7-2 


150.30 


43.84 


595-596-597 


AR-7 -2 


152.73 


42.25 


598-599-600 


AR-7 -2 


155-95 


58.68 


601-602-603 


AR-7 -2 


158. 11 


57.21 


604-605-606 


AR-7-2 


160.24 


55.74 


607-608-609 


AR-7-2 


163.84 


71.68 


610-611-612 


AR-7 -2 


165.77 


70.51 


613-614-615 


AR-7 -2 


167.72 


69. 19 


616-617-618 


AR-7-2 


169.23 


78.83 


619-620-621 


AR-7 -2 


170.60 


78.83 


622-623-624 


AR-7-2 


173.00 


7 &. 83 


625-626-627 


AR-7 -2 


168.56 


75.10 



NOTES 

1. All strain gage rosettes were mounted at the Naval 
Postgraduate School and oriented with one leg perpendicular 
to the CIB, the 45 degree gage pointing outboard and aft, 
and the third leg parallel to the CIB. 

2. All gages were manufactured by the Baldwin -Lima- 
Hamilton Corporation, taken from lot B-31 and rated at 
120.5 +0.5 ohms with a gage factor of 1.97 + 2 per cent. 

3. The high numbered rosette leg identifies the gage 
and is the leg parallel to the Load Reference Axis. 
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APPENDIX D 



STRAIN GAGE INSTRUMENTATION 

The method of taking strain readings employed common 
Wheatstone Bridge circuitry as is shown in Figure Dl. The 
experimental test setup used the Budd/Datran Digital Strain 
Indicator (Model TC22) . The Wheatstone Bridge and all 
associated circuits were incorporated within the indicator. 
The indicator interpreted the bridge unbalance and gave a 
visual display of the strain directly in units of micro- 
inches per inch. 

When the test runs were made, each of the twenty gages 
connected through the switching and balancing unit was 
zeroed on the digital strain indicator in the conventional 
manner prior to load application. At any time during the 
test run the calibration could be checked by merely pushing 
a button on the counter. This gave a continual check and 
ensured the calibration of the counter at all times during 
the running of the test. 
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FIG. Dl 

SCHEMATIC DIAGRAM 
STRAIN GAGE INSTRUMENTATION 



B 





R, = R 4 = 120.5 ohms 
R 2 = R 3 = 120.5 ohms 
Resistance tolerance: 5% 
Linearity tolerance: 0.5% 
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APPENDIX E 



EXPERIMENTAL PRINCIPAL STRESSES AND 



AXIS ORIENTATION 



The magnitudes of the principal stresses, maximum shear 
stresses and orientations were calculated for each rosette 
on both sides of the wing panel. The experimental strain 
gage readings are tabulated in Table El. Calculations were 
simplified through the use of the IBM 360 Digital Computer 
utilizing a FORTRAN IV program named "RECROSE" shown in 
Tables E2 and E3. The RECROSE results are listed in 
Table E4. 

The program was designed around the standard relation- 
ship between rectangular rosette readings and principal 
stresses found in any standard text (e.g., Wang in Refer- 
ence 7) on the subject. 



3 
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Figure El 
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<r ~- ■ T [^ - 




z£*~ fe i £*) 
£-\ * *=■ % 



(j ) p is the angle between leg number one and the maximum 

normal stress axis. A positive value indicates an angle in 

the direction of € 0 . and £_ are the strain 

^ 1 2 -i 

gage readings on their respective legs shown in Figure El. 
The number of the rosette leg parallel to the Load Reference 
Axis (CIB) is the identifying number for the rosette. 

The input to program RECROSE (Table El) was designed to 
accommodate the above notation. Only the rosette identify- 
ing number is listed. This is also the gage number for 
which the first column of strains is listed. The second 
and third columns are then the values of the strain on the 
next two consecutive numbered gages comprising the rosette. 

This program was designed for 7079-T6 Aluminum with a 
Modulus of Elasticity of 10.3 x 10^ pounds per square inch 
and a Poisson's ratio of 0.32. 

The results of RECROSE are listed in Table E4 and are 
identified by the rosette numbers which orient the reference 
axis for 0p • 
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TABLE El 



EXPERIMENTAL STRAIN READINGS 



(Input to RECROSE Program) 
336,000 in. lb. load 



INSIDE SKIN 


Rosette 


M inj 


c M Id/ 


r M it 1/ 


Number 


c, /w. 




C-3 /iv\. 


83 


9 


51 


-12 


86 


12 


59 


-10 


89 


15 


55 


-20 


92 


13 


52 


-21 


95 


12 


83 


-3 


98 


19 


83 


-8 


101 


20 


91 


-4 


104 


18 


109 


-7 


107 


10 


110 


-6 


110 


21 


107 


-12 


113 


14 


136 


-2 


116 


17 


136 


-13 


119 


22 


130 


-16 


122 


20 


147 


-11 


125 


12 


152 


-15 


128 


17 


151 


-21 


40 


26 


149 


-19 


OUTSIDE SKIN 


582 


10 


59 


-18 


585 


12 


70 


-13 


47 7 


18 


65 


-20 


588 


26 


69 


-26 


591 


22 


87 


-20 


594 


32 


92 


-26 


597 


26 


89 


-26 


600 


20 


124 


-12 


603 


25 


123 


-16 


606 


30 


120 


-18 


609 


20 


152 


-23 


612 


29 


154 


-14 


615 


22 


146 


-6 


618 


32 


169 


-19 


621 


32 


169 


-30 


624 


33 


17 7 


-13 


627 


33 


167 


-23 
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TABLE E2 





LIST OF SYMBOLS USED IN 




PROGRAM RECROSE 



Computer 
Coded Name 


Definition 


El 


£ strain in identifying leg of rosette. 


E2 


£. 2 ' strain in diagonal leg of rosette. 


E3 


£ 2 ' strain in perpendicular leg of rosette, 


SIGMAX 


rn , maximum principal stress. 


SIGMIN 


errin' minimum principal stress. 


SIGMA 1 


stress in identifying leg of rosette. 


SIGMA 3 
TAUMAX 


0” stress in perpendicular leg of rosette 
^max' n,ax i n,uin shearing stress. 


PHIPRR 


, angle from £ . to <j~ axis. 

1 ^ max 


N 


The number of input rosettes. 


E 


Modulus of Elasticity. 


U 


J/. , Poisson's ratio. 
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TABLE E3 

FORTRAN IV FROGRAM 
"RECROSE" 



//PECBOSE JCB 2802 . PARCELLS .PSGLEVEL 0 
// EXEC FCRTCLG 
//FCRT.SYSIN DO • 

RECTANGULAR STRAIN RCSETTE DATA REDUCTION FCR TFE CETERM INAT ICN 
OF PRINCIPAL AXES. PRINCIPAL STRESSES, HAXIPUP SHEARING STRESSES, 
ANC STRESSES ALCNG TFE RCSETTE AXES 1 AND 3. 



8 ft! 



READ 
L 0 
I 0 

1 I I 1 
REAC 



5,100 N,E,U 



5,101 I GAGE ,E 1 ,E2 t E 3 



El E1*1.0E-C6 
E 2 E2*1.0E-06 
E3 E3*l.QE-06 
CORRECT ICN FOR 
El E 1 * 1 .0 3 
E2 E2*1.03 
1.03 



A GAGE FACTCR 



1.95 



E3 E3* 1 
G 1.0-U 
H 1.0 U 
A El E 3 
B E1-E3 
C 2.0*E2 -A 
R SORT 0**2 C* *2 
TAUPAX E/ 2 • 0 *H *R 
SIGPAX E/2.0 * A/G R/H 

SIGPIN E/2.0 • A/G - R/H 

PHIPRR 0.5* AT AN2 C,8 
IF B.GE.O.C GC TC 2 
PHIPRR PHIPRR 1.5707563 

2 PH I PRC PHIPRR/. 1745329E-C1 

D TAUPAX*CGS 2.0*PHIPRR 
SIGPA1 SIGMAX SIGPIN /2.0 0 

SIGPA3 SIGMAX SIGPIN /2.0 -0 

TAUO -TAUPA X*S IN 2.C*PHIPRR 

IF l.GT.l GO TO 3 
WRITE 6,102 

3 WRITE 6,103 I GAG E , E 1 , S I GPAX , E2 , S IGP I N , E3 , T ALPA X , S IGMA 1 , PH I PRR , 
9SIGPA3,Ph I PRO, TAUO 

L l 11 

IF L.LT.66 GC TO 4 

WRJTE 6,102 

4 IF I.LT.N GC TO 1 
IOC FORPAT I 8 , E 17 • 7 , F 10 . 5 

101 FORPAT 18.3FI2.5 

102 FORPAT 1 H 1 

103 FORPAT T6. ROSETTE NUPBFR .I8.//,T16, El .F14.7, IN/IN ,T45, 

3 SIGPAX • H 14. 7 , PSI , / , T 16 , £2 ,F14.7. IN/IN .T45, SIGPIN , 

3F14.7, PSI a / , T 1 6 , E 3 ,F14.7, IN/IN ,t45, TAUPAX 

3 PSI ,//,Tll, SIGPA 1 ,F14.7. 

3S • / , T 1 1 , SIGMA 3 ,F14.7, PSI 

3T 1 3 , TAU C , F 14 • 7 , PSI ,// 

RETURN 
END 



, t F 1 4 . 7 , 

PSI , T45 , PHI PR .F14.7. R AC I AN 
, T45, PHI PR , F14 .7 , DEGREE S ,/, 



1 

COO 3 
C 004 
0^05 
0 006 
0007 

838? 

0 3 1C 
0:il 
OC 1 2 

00 1 3 
CO 14 

M 2 

1GFC 
2GFC 
3GFC 
4GFC 
CO 17 
0010 
OC 1 5 
C02G 
C 02 1 
0022 
0 323 
0.'24 
0" 2 5 
0326 
0027 
C'-28 
0;>29 
O ' 30 
CC 3 1 
C032 
0 33 3 
0 J34 
CO 35 

0036 

0037 
C03e 
CO 35 
034C 
C 04 1 
0042 
Q 043 
0 144 
0 04 5 
0046 
0 147 
0 8 
0C45 
CC5C 
0051 

8:11 
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TABLE E4 



EXPERIMENTAL RESULTS 
(Output of Program RECROSE) 



Ros- 

ette 


SIGMA 1 


SIGMA 3 


SIGMAX 


SIGMIN 


TAUMAX 


PHIPRR 


83 


60.987 


-107.792 


406.902 


-453.706 


430.304 


39.345 


86 


104.009 


-7 2.807 


490.063 


-458.860 


474.462 


39.630 


89 


101.646 


-179.653 


444.059 


-522.066 


483.063 


36.536 


92 


74.225 


-199.037 


407.954 


-532.765 


470.359 


36.556 


95 


130.485 


9.928 


703.993 


-563.579 


633.786 


42.271 


98 


194. 309 


-22.693 


718.064 


-546.448 


632.256 


40.059 


101 


221.257 


28.366 


798.828 


-549.205 


674.016 


40.886 


104 


186.273 


-14.657 


923.694 


-752.078 


837.886 


41.556 


107 


95.500 


-33.094 


901.589 


-839.184 


870.386 


42.881 


110 


202.819 


-62.407 


904.616 


-7 64. 203 


834.409 


40.427 


113 


157.906 


29.310 


1140.408 


-953.191 


1046.799 


43.239 


116 


151.760 


-89.355 


1114.900 


-1052.495 


1083.697 


41.806 


119 


199.510 


-105.901 


1078.876 


-985.268 


1032.072 


40.745 


122 


194.782 


-54.369 


1222.249 


-1081.836 


1152.042 


41.896 


125 


85.100 


-131.904 


1215.054 


-1261.858 


1238.457 


42.487 


128 


121.503 


-183.908 


1207.919 


-1270.324 


1239.123 


41.461 


40 


235.441 


-126.231 


1237.903 


-1128. 693 


1183.298 


40.605 


582 


50.114 


-174.925 


456.284 


-581.095 


518.689 


38.736 


585 


92.664 


-108.265 


567.653 


-583.255 


575.454 


39.973 


477 


137.104 


-168.307 


536.391 


-567.594 


551.992 


36.970 


588 


208.965 


-208. 965 


592.624 


-592.624 


592. 624 


34.646 


591 


184. 381 


-153.178 


727. 102 


-695.899 


711.500 


38.139 


594 


279.881 


-186.272 


799.123 


-705.514 


752.318 


35.976 


597 


208.965 


-208.965 


745.201 


-745.201 


745.202 


36.857 


600 


191.000 


-66.189 


1035.393 


-910.581 


972.987 


41.203 


603 


234.968 


-94.555 


1036.749 


-896.336 


966.543 


40.093 


606 


286.500 


-99.282 


1029.923 


-842.706 


936.314 
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